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Photoproduction of heavy quarks in uhraperipheral colhsions can help elucidate important fea- 
tures of the physics of heavy quarks in Quantum Chromodynamics (QCD). Due to the dependence 
on parton distributions it can also potentially offer some constraining ability in the determina- 
tion of nuclear parton distributions. In the present study we consider next-to-leading order (NLO) 
photoproduction of heavy quarks in ultraperipheral proton-proton (pp), proton- nucleus (pA), and 
nucleus-nucleus (AA) collisions at the CERN Large Hadron Collider (LIfC). Total cross sections 
and rapidity distributions are considered and the influence of nuclear modifications of parton distri- 
butions on these quantities are explored for pA and A A collisions. We find that photoproduction of 
heavy quarks in PbPb collisions exhibit significant sensitivity to nuclear effects, and in conjunction 
with photoproduction in pPb collisions, aff'ords good constraining potential for gluon shadowing 
determination. 
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I. INTRODUCTION 

The study of ultraperipheral relativistic heavy ion col- 
lisions is useful in exploring several aspects of particle 
and nuclear physics and is an important part of current 
experimental efforts at the LHC. Consequently ultrape- 
ripheral collisions have been extensively discussed in the 
literature (for a small sample of references, see e.g. [J- 
fisj ) for different collision systems and various photon- 
nucleus and photon-photon processes. 

In [l^ - [l6j we have considered photoproduction of cc 
and bb in ultraperipheral pPb (^s„„ — 5 TeV and 
(^/i;;^ = 8.8 TeV) and PbPb {./s^ = 2.76 TeV and 



8.8 ievj and fbtb (^s, 
iy/^NN — 5-5 TeV) collisions at the LHC. A major short- 



coming was that the results were calculated to leading or- 
der (LO) in perturbative QCD. In this article we remedy 
this shortcoming and present results at next-to-leading 
order (NLO) accuracy. We also consider, as an addition 
to the previous treatments, photoproduction of heavy 
quarks in ultraperipheral pp collisions at y/s^ = 7 TeV 
and ^s„„ =14 TeV. As in the earlier studies an im- 
portant goal is to investigate quantitatively the influence 
of nuclear modifications of parton distributions on ob- 
servables such as cross sections and rapidity distributions 
and also the extent of applicability of these observables 
in helping to constrain components of nuclear parton dis- 
tributions. 

The paper is organized as follows: in Sec.|TT]we discuss 
photoproduction of heavy quarks and the basic ingredi- 
ents required in our calculations, namely photon flux, 
photohadron cross sections, and input nuclear and pho- 
ton parton distributions. In Sec. IIIII we present the re- 
sults of our calculations. Our conclusion is contained in 
SecHVl 



II. PHOTON-HADRON INTERACTIONS IN 
ULTRAPERIPHERAL COLLISIONS 

In general the total cross section for the photopro- 
duction of a pair of heavy quarks (cc or bb) in ultra- 
peripheral collisions can be expressed as a convolution 
of the incident photon flux, dNj/dk, with the energy- 
dependent photo-cross section a''^^'^'^ (k), with k de- 
noting the photon energy and H a proton or nucleus. 
Thus for pp collisions we can write 
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with dNP / dk the photon flux from a proton in the pp col- 
lision system and a^^^^^^ (k) the energy-dependent cross 
section for the photoproduction of a QQ pair off a proton. 
The factor of 2 accounts for the source/target symmetry 
present in symmetric collisions, since each proton can 
act simultaneously as a source and target of photons. A 
similar expression holds for AA collisions with dN^/dk 
replaced with dN^ /dk, the flux from one of the partici- 
pating nuclei and a^^^'^'^ (k) replaced with a'^^^'^'^ {k), 
the cross section for photoproduction of a QQ pair off 
a nucleus. In the case of pA collisions the nucleus acts 
preferentially as the source and the proton as the target, 
leading predominantly to 7p processes. But there is still 
a non-negligible contribution from 7A processes in which 
the proton acts as the source of photons and the nucleus 
as the target. Thus expressions for both types of fluxes 
are required for pA collisions. We can thus write 
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with dN^ /dk and dN'P/dk the fluxes from the nucleus 
and proton in the pA colhsion system and a^P^'^'^ (k) 
and a'^^^^'^{k) the photonucleon and photonuclear 
cross sections respectively. 

Manifestations of nuclear effects can be more trans- 
parently observed when rapidity distributions are con- 
sidered. Using the relation da/dy = kda/dk the pair 
rapidity distribution can be expressed as 
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and scales directly with the photon flux dN^/dk. Thus 
for pA collisions, with the convention that the proton is 
incident from the right and the nucleus from the left, the 
total rapidity distribution is 
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where fc; (A;; cx e~^) and fc^ (kr oc e^) simply denote pho- 
tons from the nucleus and proton respectively. The first 
term on the right-hand side (7p distribution) peaks at 
positive rapidities while the second term (7A distribu- 
tion) peaks at negative rapidities. Since both the fluxes 
and cross sections are different, the total distribution is 
manifestly asymmetric, and the 7p term dominates due 
to the much larger nuclear flux dN^ /dk. 

The total rapidity distribution for AA collisions can 
likewise be written as 
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with ki and kr simply denoting photons from the nucleus 
incident from the left and right respectively. Here the 
cross sections and the left/right fluxes are identical; thus 
the flrst term on the right-hand side of Eq. ([5]) yields a 
distribution which is the mirror image of the second term. 
Consequently the total distribution is symmetric about 
midrapidity {y = 0). Similar expression and symmetry 
attribute hold true for pp collisions. 
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which drops off exponentially for C > 1, above a cutoff 
energy determined essentially by the size of the nucleus, 
Sc«to// ~7MeV/6 (fm). 

The photon flux from a proton is usually estimated 
using the dipolc formula for the electric form factor QJ^ : 
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where £1 = 1 + [0.71 GeVVQ^^j 
mentum transferred Q^jj^ = fc^/[7^(l 

In the case of proton-nucleus (pA) collisions, the flux 
from the proton is given by Eq. ([7]) while the flux due to 
the nucleus (of charge Z) can be evaluated analytically 
and is given by [2] , 



dN^jk) 
dk 



2Z^a 
Trk 



■pA\ 



(C: 



pA\2 



2(/-pAy^ 



(8) 



with reduced adiabaticity parameter, C,^ , given by 
— k{Rp + Ra)/i and Rp the effective radius of the 
proton. 

For nucleus- nucleus (AA) collisions an analytic approx- 
imation similar to Eq. ([5]) can be derived. A more accu- 
rate expression for the flux can be obtained by integrat- 
ing d'^N^{k,h)/dkd^b over impact parameters with the 
constraint of no hadronic interactions and accounting for 
the photon polarization. This yields the total photon flux 
dN^{k)/dk given by 
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with Ra the radius of the nucleus. 

B. Photoproduction cross sections 



(9) 



A. Photon flux 

Let us now consider the various photon fluxes oc- 
curring in the expressions for total cross sections and 
rapidity distributions. For a given impact parameter 
b, the flux of virtual photons with photon energy fc, 
d^ Nj{k,h) / dkd'^b, is strongly dependent on the Lorentz 
factor 7. The photon flux also depends strongly on the 



In general photon interactions with hadrons and nuclei 
can be classifled as direct or resolved. In direct interac- 
tions the photon behave as a point-like particle while in 
resolved interactions the incident photon first fiuctuates 
into a quark- antiquark state (or an even more complex 
partonic configuration consisting of quarks and gluons) 
which then subsequently interacts hadronically with the 
hadron or nuclear target. The cross section for the pho- 
toproduction of a pair of heavy quarks, crT^^^-QQ jg 
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therefore a sum of both the direct and resolved contribu- 
tions: 



<^'- - ik) = -ntr ik)+-:fsl'i? (k). (10) 



Here H stands for a proton or a nucleus {H = p, A) 
and the total photoproduction cross section is obtained 
by convolution of the energy-dependent cross section 
^-yH^QQ g^^^ i-j^g equivalent photon flux, dN^{k)/dk, 
as given in Eq. ([T]) for symmetric pp (AA) collisions and 
Eq. ([2]) for asymmetric pA collisions. 

Let us consider both contributions in some detail. In 
view of the high energies involved perturbative QCD 
is applicable, and both the direct and resolved contri- 
butions are expressible as convolution of the cross sec- 
tions for the relevant partonic subprocesses and the cor- 
responding nucleon/nuclear parton distributions. Thus 
the direct contribution can be written as 



direct 



x/f(x„,Q2)e(c.). (11) 



Here ^'^"^'^'^(xqS, Q^) is the parton-level cross section 
for the photoproduction of a heavy quark pair from the 
interaction of a photon 7 and a parton a with momentum 
fraction Xa- The renormalization scale, /^r, has been set 
equal to the factorization scale, /iy i.e. jjLr = jJLf = Q 
and {xa, Q^) is the parton distribution of a in eval- 
uated at Xa and Q^. In addition s — denotes the 
square of the center-of-mass energy of the photon-nucleus 
(or photon- nucleon) system, s = W^^ — XaS that of the 
photon-parton system, and Ca = s — 4r7iQ with toq the 
mass of the heavy quark (charm or bottom). The func- 
tion 0(Ca) enforces a minimum ("threshold") value of Xa, 
xj^*", on the integral given by x™™ = Am^/s. Note that 
the summation over q involves only the light flavors, i.e. 
q = u,d, s. 

The resolved contribution is similar in structure to 
hadroproduction of heavy quarks and can be written as 
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where a"''''^'^'^ {xaXhS, Q^) is the partonic ab — > QQ cross 
section, f2{xa,Q'^) [f^ {x2,Q'^)) is the distribution of 
parton a (b) with momentum fraction Xa (xb) in a photon 
(H) respectively, s — XaXbS and Cab — s ~ Avoq. Similar 
to the direct contribution, the summation over q involves 
only light quark flavors, i.e. q — u,d, s. 

It is thus clear that the two major ingredients needed 
for the determination of the hadronic photoproduction 
cross sections are the partonic (parton-level) cross sec- 
tions and parton distributions in nucleon/nuclei and pho- 
tons. Below we briefly discuss them. 



1. Partonic cross sections 

Photon-gluon fusion is the only relevant subprocess for 
the direct contribution at leading order (LO), and the ex- 
pression for the cross section can be found in (T8l - [20| . In 
the case of LO resolved contribution, only the gg and qq 
channels are relevant, and the corresponding cross sec- 
tions can be found in [2l| - |23j . We have used the next- 
to-leading order (NLO) results presented in [23, HH] for 
both direct and resolved contributions. Small- a; effects 
(see [2^) which could potentially be important for charm 
photoproduction are not addressed in the present study. 
These effects should be somewhat lessened in the case of 
the resolved contribution due to the inherent small-x cut- 
off {x > 0.001) present in the photon parton distributions 
used in this work. 



2. Parton distributions tn nuclet and photons 

We now turn to the consideration of the parton distri- 
butions (PDs) relevant to the processes considered in the 
present study. While the direct contribution to the pho- 
toproduction of heavy quarks requires the distributions of 
light quarks/antiquarks and gluons in protons and nuclei, 
the resolved contribution, in addition, requires these dis- 
tributions in photons also. Thus in addition to the usual 
requirement of nucleon and nuclear parton distributions 
(nPDs), there is also the need for the relatively poorly 
known parton distributions in photons (7PDS), thereby 
increasing the level of the theoretical uncertainties in the 
calculation of photoproduction of heavy quarks. 

Let us first discuss nuclear parton distributions. It is 
rather well-known that the distributions of partons (i.e. 
quarks and gluons) in nuclei are quite different from the 
distributions in free nucleons, that is, they are "modi- 
fied" by the complex, many-body effects in the nucleus. 
These nuclear effects are usually parametrized in terms 
of "nuclear modifications" i?^(x, Q^) which in general 
depend on the parton specie (a), the nucleus {A), mo- 
mentum fraction x and scale . The nuclear effects can 
be categorized based on different intervals in x. At small 
values oi x [x < 0.04), we have the phenomenon gener- 
ally referred to as shadowing. This is a depletion, in the 
sense that in this interval, the distribution of a parton a 
in the nucleus is smaller compared to the corresponding 
distribution in a free proton, i.e. < 1. Antishadow- 
ing, which is an enhancement {R^ > 1), occurs in the 
range 0.04 < a; < 0.3. Another depletion, the classic 
EMC effect 27], is present in the interval 0.3 < a; < 0.8, 
while for x > 0.8, the Fermi motion region, we have an- 
other enhancement. It is important to note that although 
both shadowing and the EMC effect (antishadowing and 
Fermi motion) correspond to depletion (enhancement), 
the physical principles and mechanisms governing these 
phenomena are quite different. Further details can be 
found in [21111 With the knowledge of R^{x,Q^), nu- 
clear parton distributions can be expressed as a convo- 
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lution of free nucleon parton distributions and nuclear 
modifications, i.e. f^{x,Q^) = faix,Q^) (8) R^{x,Q^). 

While the determination of quark and antiquark distri- 
butions in nucleons and nuclei is in general a nontrivial 
task, that of gluons is even more problematic. Gluons 
are electrically neutral, and thus their distributions can- 
not be extracted directly from Deeply Inelastic Scattering 
(DIS) and Drell-Yan (DY) processes which account for 
the major part of the data used in global fits. Their dis- 
tributions are in general inferred from sum rules and the 
evolution of sea quarks distributions. The situation is 
even worse in the nuclear case: the available data is much 
less than for nucleons, and there is the added complica- 
tion of a mass dependence. It is therefore not unusual 
for nuclear gluon distributions from different global fits 
to differ significantly, especially in the magnitude of the 
various nuclear effects (shadowing, antishadowing, etc). 
This is especially obvious at low (i.e. around their 
initial starting scales) since evolution to high tends 
to lessen the differences. Earlier global analyses [32I - 
[ssj relied heavily on fixed-target nuclear deep-inelastic 
scattering (DIS) and Droll- Yan (DY) lepton-pair produc- 
tion data. Incorporation of additional data from inclu- 
sive hadron production in deuteron-gold collisions and 
neutrino- iron processes has been implemented in 36 4l| . 
The approach in [i^, utilizes the Gribov picture of 
shadowing and data from diffractive processes to gen- 
erate nuclear modifications. Despite all these advances 
the nuclear gluon distribution is still currently the least 
constrained aspect of nuclear parton distributions, as sig- 
nificant uncertainties still persist at both small and large 

X. 

Three recent nucleon/nuclear parton distributions are 
utilized in the present study. For the proton we use the 
Martin-Stirling-Thorne- Watts (MSTW08) parton distri- 
butions 44] which are available up to next-to-next-to- 
leading order (NNLO). In the nuclear case we use two nu- 
clear modification sets: the NLO set by Eskola, Paukun- 
nen, and Salgado (EPS09) [13 and the NLO set by 
Frankfurt, Guzey, and Strikman (FGSIO) Q. For 
EPS09 we use the central fit while for FGSIO we use 
the strong gluon shadowing model (FGSIO-H). The dis- 
tributions from MSTW08 serve two purposes: as the free 
nucleon distributions used in conjunction with nuclear 
modifications, and also as a "special" nuclear distribu- 
tion in the absence of nuclear effects. The latter case is 
particularly useful for highlighting the influence of the 
various nuclear effects on observables. 

It is instructive to compare the characteristics of the 
light partons and gluon distributions from both EPS09 
and FGSIO sets based on the strength of their nuclear 
modifications. In Fig. [T] we show the nuclear modifica- 
tions for light quarks and gluons in Pb, Rjb{x, Q^) from 
EPS09 and FGSIO at the factorization scale = Ami. 
For valence up (m„) and down (d^) quarks, both distri- 
butions exhibit identical antishadowing, EMC, and Fermi 
motion characteristics. Both and are less shadowed 
in FGSIO than in EPS09 for the entire shadowing interval 
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FIG. 1. (Color online) Nuclear modifications of light quarks 
and gluons in Pb, Rf{x, = Ami GeV^), from EPS09 (solid 
line) and FGSIOH (dashed line): (a) valence up quark u^, (b) 
up sea quark Ua (= u), (c) valence down quark d^, (d) down 
sea quark ds (= d), (e) strange quark s (= s), and (e) gluon 

g- 



(w 0.9 versus 0.8 at very small x). In the case of u and 
d, the EMC and Fermi motion effects which are present 
in EPS09 are virtually absent in FGSIO. The shadow- 
ing components are not very different, with FGSIO being 
slightly more shadowed at very small x, and EPS09 more 
shadowed for x > 10^"*. In FGSIO strange quarks (an- 
tiquarks) are more shadowed but there is an absence of 
EMC and Fermi motion characteristics which are present 
in EPS09. By far the most impressive difference occurs in 
the shadowing of the gluon sector: the gluon shadowing 
in EPS09 is almost constant at small x at w 0.8 while in 
FGSIO the shadowing increases with decreasing x, reach- 
ing a value of about 0.49 at x « lO'^. While EPS09 
exhibits both EMC and Fermi motion effects these are 
absent in FGSIO. Also the antishadowing component is 
smaller in FGSIO than in EPS09. Thus in the shadowing 
regime it is reasonable to expect that much of the dif- 
ferences in the predictions from EPS09 and FGSIO stem 
from their radically different gluon characteristics. 

Parton distributions in photons (7PDS) are derived 
from experimentally determined photon structure func- 
tion F2 {x^Q'^), in conjunction with appreciable theoret- 
ical inputs. These inputs, which are necessary in im- 
plementing the parametrization of photon parton distri- 
butions from the structure function, account in part for 
some of the observable differences in the available pho- 
ton parton distribution sets. Another source of differ- 
ences is in the choice and scope of experimental data from 
which F2 is extracted. At present there is an apprecia- 
ble number of photon parton distribution sets available, 
both at leading and next-to- leading orders (45l - l6]| . It 
should be noted that unlike in the case of a nucleon, there 
are no valence quarks present in the photon; therefore 
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antiquark distributions are the same as quark distribu- 
tions. Furthermore, there are no sum rules governing the 
gluon content; thus the gluon distribution is almost to- 
tally unconstrained. The gluon distribution contributes 
to F2{x,Q^) majorly through the 7*3 — > qq channel, 
which has significant numerical support only at small x. 

Two recent NLO photon parton distributions are the 
set by Cornet, Jankowski, and Krawczyk (CJK04) [g^l 
and that by Aurenche, Fontannaz, and Guillet (AFG04) 
[6l| . The C JK04 set uses the DIS renormalization scheme 
while the MSbar scheme is employed in AFG04. We have 
used the AFG04 set (AFG04_BF) in the present study 
for consistency with the renormalization scheme used in 
[13, HE] for the NLO cross sections used in this work. In 
Fig.[2]we show the parton distributions from the AFG04 
NLO set for the light (u,d,s) quarks and gluon. The most 
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FIG. 2. (Color online) Parton distributions in the photon 
for (a) up, (b) down, (c) strange quarks, and (d) gluons at 
g2 = /^rric GeV^ from the AFG04 distribution set. For better 
visuality the gluon distribution (d) has been scaled down by 
a factor of 10. 

remarkable characteristic is the dramatic rise of g'^ with 
decreasing x. 

III. RESULTS 

We now discuss the results of our calculations for the 
inclusive photoproduction of heavy quarks (cc and bb) in 
ultraperipheral proton-proton (pp), proton-lead (pPB), 
and lead-lead (PbPb) collisions at the LHC. The calcu- 
lations are carried out at two center-of-mass (cms) en- 
ergies for each collision system: E'f {y/s^ — 7 TeV) 
and i^W^ = 14 TeV) for pp, Ef = 5 

TeV) and Ef'' (^s^ = 8.8 TeV) for pPb, and E^''^'' 
(ys^ = 2.76 TeV) and E^^'' (^/5^ = 5.5 TeV) for 
PbPb collisions respectively. We take rric = 1.4 GeV and 
nib = 4.75 GeV for consistency with the MSTW08 par- 
ton distributions. The strong coupling constant at scale 



Q^, as{Q'^), is evaluated to two loops using the evolution 
code contained in the MSTW08 package. For cc we set 

A. Heavy quarks in pp collisions 

In Table [T] we show the cross sections for the direct and 
resolved components, as well as their sum, for the pho- 
toproduction of cc and bb in ultraperipheral pp collisions 
at the two specified cms energies. 

TABLE I. Direct, resolved, and total cross sections for pho- 
toproduction of cc and bb in ultraperipheral pp collisions at 
the LHC. All cross sections are in nanobarns (nb). 





(TeV) 


Direct 


Resolved 


Total 


cc 


7 


2050.9 


677.1 


2728.0 




14 


2837.6 


1124.7 


3962.3 


bb 


7 


31.2 


20.1 


51.3 




14 


52.8 


44.0 


96.8 



Some interesting features can be deduced from Table HI 
Thus at y/s^j^ — 7 TeV the resolved contribution is about 
24.8% of the total cross section for cc and 39.1% for bb, 
while at ^s„„ = 14 TeV the resolved component ac- 
counts for 28.4% and 45.5% respectively for cc and bb 
production. It is thus evident that the resolved contribu- 
tion, being a significant fraction of the total cross section, 
is important for cc, and more especially so for bb pro- 
duction, with this importance increasing with increasing 
energy. 

These features are apparent in the rapidity distribu- 
tions for cc as shown in Fig. [3l and for bb as displayed 
in Fig. m These distributions are manifestly symmet- 
ric about midrapidity {y = 0) due to the source/target 
symmetry present in pp collisions. 

It is instructive to explore how the cross sections scale 
with energy. To this end we consider the ratio of the to- 
tal cross section at ^s„„ = 14 TeV to that at ^s„„ = 7 
TeV for both cc (denoted by Rcs) and bb (denoted by R^i) 
production. From Table U we can see that i?cc — 1.452 
and Rfji = 1.887 respectively. Thus for cc the total 
cross section at ^/s^ = 14 TeV is approximately 1.5 
times larger than at ^/s^ = 7 TeV, and approximately 
1.9 times larger in the case of bb production. This im- 
plies that the cross section for bb production grows more 
rapidly with increasing energy. Componentwise the ratio 
of direct cross sections is « 1.4 for cc and « 1.7 for bb 
while for resolved it is « 1.7 and « 2.2 respectively. 

B. Heavy quarks in pPb collisions 

Photoproduction in ultraperipheral pPb collisions arise 
from two distinct processes: 7p interactions in which the 
Pb nucleus is the source of the photon and the proton 
acts as the target, and 7Pb interactions in which the 
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FIG. 3. (Color online) Rapidity distributions of cc photo- 
production in ultraperipheral pp collisions at the LHC. Solid 
lines in (a), (c), and (e) depict the direct, resolved, and total 
distributions at ys„„ = 7 TeV while solid lines in (b), (d), 
and (f) are the equivalent distributions at ys„„ = 14 TeV. 
In the two upper rows dashed and dash-dotted lines denote 
contributions from photons incident from the left and right 
respectively. 
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FIG. 4. (Color online) Rapidity distributions of hb photo- 
production in ultraperipheral pp collisions at the LHC. Solid 
lines in (a), (c), and (e) depict the direct, resolved, and total 
distributions at -Js^ = 7 TeV while solid lines in (b), (d), 
and (f) are the equivalent distributions at y^s^^iv ~ 1^ TeV. 
In the two upper rows dashed and dash-dotted lines denote 
contributions from photons incident from the left and right 
respectively. 



proton is the source of photons and the Pb nucleus is the 
target. Due to the large disparity in the magnitude of the 
photon fluxes from the proton and the nucleus (Pb) the 
7P contribution dominates, and a common practice is to 
neglect the photonuclear (7Pb) contribution. We retain 



the photonuclear contribution in this study and test the 
severity of the approximation made in neglecting it. 

In Table [TTl we present the cross sections for the direct 
and resolved components of both 7p and 7Pb contribu- 
tions for cc and 66 production in ultraperipheral pPb col- 
lisions at the LHC. Total cross sections for both cc and 
66 are presented in Table Hill 



TABLE II. Direct and resolved cross sections for photopro- 
duction of cc and bb in ultraperipheral pPb collisions at the 
LHC. All cross sections are in microbarns (fih). 











= 5 TeV 




= 8.8 TeV 






PDF 


Direct 


Resolved 


Direct 


Resolved 




7P 


MSTW08 


2474.6 


462.2 


3563.5 


823.8 


cc 




MSTW08 


179.6 


53.8 


238.2 


84.0 




7Pb 


EPS09 


157.9 


48.7 


207.2 


74.4 






FGSIO 


131.9 


46.5 


169.5 


70.5 




7P 


MSTW08 


21.4 


6.4 


38.0 


15.1 


bb 




MSTW08 


2.5 


1.4 


3.9 


2.7 




7Pb 


EPS09 


2.3 


1.4 


3.6 


2.6 






FGSIO 


2.1 


1.3 


3.2 


2.6 



Let us first consider cc production. At ^s„„ — 5 
TeV the 7Pb contribution to the total cross section is 
approximately 7.4% for MSTW08 (no nuclear modifi- 
cations), 6.6% for EPS09 (moderate gluon shadowing), 
and 5.7% for FGSIO (strong gluon shadowing) respec- 
tively. The resolved component accounts for 16.3% of 
total cross sections, and is independent of the choice of 
parton distribution set. Overall the sensitivity to nuclear 
modifications (dominantly shadowing) is very small: the 
no-modification (MSTW08) total cross section is reduced 
by approximately 0.9% and 1.7% by the modifications in 
EPS09 and FGSIO respectively. 

The trend is similar at y/s^ ~ 8.8 TcV. Here the 
7Pb contribution to the total cross section is approxi- 
mately 6.8% for MSTW08, 6% for EPS09, and 5.2% for 
FGSIO respectively. The resolved component yields ap- 
proximately 19.3% of total cross sections, and is also in- 
dependent of the choice of parton distribution set. Iden- 
tically the sensitivity to nuclear modifications remains 
the same as at ^s„„ — 5 TeV. 

These features are reflected in Fig. [5] where we show 
the component as well as total rapidity distributions at 
y/s^^ — 5 TeV and ^s„„ = 8.8 TeV respectively. In ac- 
cordance with our convention the 7p contribution peaks 
at positive rapidities while the 7Pb distributions peak at 
negative rapidities, and the asymmetric nature of the to- 
tal distributions is clearly exhibited. Due to the smallness 
of the 7Pb contribution, the effects of nuclear modifica- 
tions are totally negligible in the total rapidity distribu- 
tions as the distributions from the three different parton 
distribution sets practically overlap. 

We now discuss the case for 66 production. In line 
with the treatment of cc production, we determine the 
relative importance of 7Pb and resolved contributions to 
the total 66 production cross sections. At ^s„„ = 5 TeV 
the 7Pb contribution to the total cross section is approx- 
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FIG. 5. (Color online) Rapidity distributions of cc photo- 
production in ultraperipheral pPb collisions at the LHC. In 
(a), (c), and (e) we show the direct, resolved, and total dis- 
tributions at .^s„„ = 5 TeV while (b), (d), and (f) dis- 
play the equivalent distributions at ys„j^ = 8.8 TeV. In the 
first and second rows long-dashed line depicts the 7p compo- 
nent while short-dashed (MSTW08), solid (EPS09), and dash- 
dotted (FGSIO) lines correspond to 7Pb component with no 
shadowing, moderate shadowing, and strong shadowing re- 
spectively. In the third row short-dashed, solid, and dash- 
dotted hues depict total distributions using MSTW08, EPS09, 
and FGSIO parton distributions respectively. 



imately 12.1% for MSTW08, 11.6% for EPS09, and 11% 
for FGSIO respectively. It is immediately apparent that 
at this energy the 7Pb contributions are more significant 
than their counterparts in cc production, and also less 
sensitive to nuclear effects. The resolved component is 
responsible for 24.6% of total cross sections and is in- 
dependent of the choice of parton distribution set. The 
sensitivity to nuclear modifications is smaller than in cc 
production: the no-modification (MSTW08) total cross 
section is reduced by approximately 0.05% and 1.2% re- 
spectively by the modifications in EPS09 and FGSIO par- 
ton distribution sets. 



TABLE III. Total cross sections for photoproduction of cc 



and bb in ultraperipheral pPb collisions at Ef 



pPb 



TeV) and Ef"" 
microbarns (/ib). 



= 5 



8.8 TeV). All cross sections are in 



PDF 


cc 


bb 




Er' R,i 


MSTW08 

EPS09 

FGSIO 


3170.3 4709.5 1.486 

3143.4 4668.9 1.485 
3115.3 4627.3 1.485 


31.6 59.8 1.894 
31.4 59.4 1.891 
31.2 59.0 1.892 



Similar features are present at ^/s^ = 8.8 TeV. The 
7Pb contribution to the total cross section here is ap- 
proximately 11.1% for MSTW08, 10.4% for EPS09, and 
9.8% for FGSIO respectively. Again these values are rel- 



atively more significant than those obtained in cc pro- 
duction. The resolved component accounts for approx- 
imately 30% of total cross sections, and is again inde- 
pendent of the choice of parton distribution set. The 
sensitivity to nuclear modifications are similar to those 
obtained at Ts;;^ = 5 TeV: 0.07% for EPS09 and 1.4% 
for FGSIO respectively. 
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FIG. 6. (Color online) Rapidity distributions of 6& photo- 
production in ultraperipheral pPb collisions at the LHC. In 
(a), (c), and (e) we show the direct, resolved, and total dis- 
tributions at ^s^^ = 5 TeV while (b), (d), and (f) dis- 
play the equivalent distributions at ys„„ = 8.8 TeV. In the 
first and second rows long-dashed line depicts the 7p compo- 
nent while short-dashed (MSTW08), solid (EPS09), and dash- 
dotted (FGSIO) lines correspond to 7Pb component with no 
shadowing, moderate shadowing, and strong shadowing re- 
spectively. In the third row short-dashed, solid, and dash- 
dotted lines depict total distributions using MSTW08, EPS09, 
and FGSIO parton distributions respectively. 

In Fig.|n]we show the component as well as total rapid- 
ity distributions for hh production at ^Js^, — 5 TeV and 



8.8 TeV respectively. These distributions faith- 



fully manifest the features described above. As in cc pro- 
duction the 7p contribution peaks at positive rapidities 
while the 7Pb distributions peak at negative rapidities, 
the total distributions being thus asymmetric due to the 
dominance of the 7p contribution. The effects of nuclear 
modifications on total rapidity distributions are totally 
negligible, as the distributions from the three different 
parton distribution sets practically overlap. This is due 
to the fact that even though the 7Pb contributions are 
relatively larger, they are also almost totally insensitive 
to nuclear effects. 

Let us now discuss the magnitude of the change in cross 
section with energy. In Table IIIII we present the ratio of 
the total cross section at y/s^ — 8.8 TeV to that at 
^/s'^ — 5 TeV for cc, Rcs, and for bb, R^^. As can be 
seen from Table ITTTl j?^^ « 1.5 and i?^;, « 1.9 for all three 
parton distribution sets. These values are practically the 
same as obtained for pp collisions. 
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C. Heavy quarks in PbPb collisions 

We now discuss production of cc and bb in ultraperiph- 
eral PbPb collisions at energies y/s^ =2.76 TeV and 
= 5.5 TeV. In Table ITVl we show the cross sections 
for both direct and resolved components while total cross 
sections are presented in Table |V] 

TABLE IV. Direct and resolved cross sections for photopro- 
duction of cc and 6b in ultraperipheral PbPb collisions at the 
LHC. All cross sections are in millibarns (mb). 









= 2.76 TcV 




= 5.5 TeV 




PDF 


Direct 


Resolved 


Direct 


Resolved 




MSTW08 


552.6 


69.5 


926.2 


161.5 


cc 


EPS09 


517.6 


69.6 


845.9 


156.9 




FGSIO 


475.7 


67.6 


755.4 


152.0 




MSTW08 


3.24 


0.60 


7.52 


2.06 


bb 


EPS09 


3.30 


0.64 


7,40 


2.14 




FGSIO 


3.18 


0.62 


7.10 


2.10 



Due to both participants being nuclei, nuclear effects 
are more pronounced in PbPb collisions than in pPb col- 
lisions. Thus the resolved components, which in pPb 
collisions are quite independent of the choice of parton 
distribution set, and so by implication insensitive to nu- 
clear effects, now acquire some sensitivity in PbPb colli- 
sions. Total cross sections and rapidity distributions also 
exhibit enhanced sensitivity. 

As usual let us first consider cc production. At 
y/s^ = 2.76 TeV the resolved components are 11.2%, 
11.9%, and 12.5% respectively for MSTW08, EPS09, and 
FGSIO distribution sets. In addition, the MSTW08 total 
cross section is reduced by 5.6% and 12.7% respectively 
by the nuclear modifications in EPS09 and FGSIO sets. 

The same trend is observed at y/s^ =5.5 TeV. Here 
the resolved components are 14.9% (MSTW08), 15.6% 
(EPS09), and 16.8% (FGSIO) respectively. The influence 
of shadowing is more appreciable; the no-modification 
MSTWQ8 total cross section is reduced by 7.8% (EPS09) 
and 16.6% (FGSIO) respectively. 

These features are manifested in the rapidity distribu- 
tions shown in Fig. [T) The distributions are symmetric 
about midrapidity {y = 0) as expected for symmetric 
collision systems. Also the sensitivity to nuclear mod- 
ifications is more transparent here than in total cross 
sections. 

Shadowing is the dominant nuclear effect for —2.5 < 
y < 2.5 at = 2.76 TeV and for -3 < y < 3 at 

^s„„ = 5.5 TeV. The rapidity distributions in these 
intervals reproduce the observed trend of gluon shad- 
owing strength exhibited in Fig. [1] MSTW08 with its 
zero gluon shadowing gives the largest rapidity distri- 
bution while FGSIOH, with its strong gluon shadowing, 
gives the smallest. Due to strong flux suppression, shad- 
owing is most markedly apparent for the rapidity range 
-1-5 < y < 1-5 for ^/s;;^ = 2.76 TeV and -2 < ?/ < 2 
for ^s„„ =5.5 TeV. These intervals are therefore partic- 
ularly suited for constraining purposes and also for dis- 
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FIG. 7. (Color online) Rapidity distributions of cc photopro- 
duction in ultraperipheral PbPb collisions at the LHC. In (a), 
(c), and (e) we show the direct, resolved, and total distribu- 
tions respectively at = 2.76 TeV while (b), (d), and (f) 
denote the equivalent distributions at y^s^^iv" ~ 5.5 TeV. In 
all cases long-dashed (MSTW08), solid (EPS09), and dash- 
dotted (FGSIO) lines correspond to no shadowing, moderate 
shadowing, and strong shadowing respectively. 



criminating among different gluon shadowing scenarios. 

At y/s^ =2.76 TeV the rapidity intervals 2.5 < y < 5 
corresponds to Xmm in the antishadowing region (deep 
shadowing) for right (left) incident photons and vice 
versa for— 5<j/<— 2.5. Due to the photon flux suppres- 
sion in the deep shadowing region, the rapidity distribu- 
tions are sensitive mainly to antishadowing in addition to 
both EMC effect and Fermi motion. Since both EPS09 
and FGSIOH have substantial antishadowing, their ra- 
pidity distributions reflect this, being slightly larger than 
those from MSTW08. The discriminatory power here 
is not as appreciable as in the shadowing case though, 
due largely to the smallness of the distributions. This 
is also the case for the rapidity intervals 3 < y < 6 and 
-6 < y < -3 at = 2.76 TeV. 

TABLE V. Total cross sections for photoproduction of cc and 
bb in ultraperipheral PbPb collisions at Ei^'^'' (.ys„„ = 2.76 
TeV) and _£;|'''^'' (^/s]^ = 5.5 TeV). All cross sections are in 
millibarns (mb). 



PDF 


cc 


bb 


Hi 1^2 ticc 


^1 ^2 ^bb 


MSTW08 

EPS09 

FGSIO 


622.0 1087.7 1.75 

587.2 1002.8 1.71 

543.3 907.3 1.67 


3.9 9.6 2.49 
3.9 9.5 2.42 
3.8 9.2 2.41 



Let us now consider the corresponding case of bb pro- 
duction in ultraperipheral PbPb collisions. The sensi- 
tivity to nuclear modiflcations, while more appreciable 
here than in pPb collisions, is signiflcantly less than for 
the equivalent cc production. At ^/s^ = 2.76 TeV the 
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resolved components are 15.7%, 16.2%, and 16.3% re- 
spectively for MSTW08, EPS09, and FGSIO distribution 
sets. Thus to a good approximation the resolved compo- 
nent can be taken as 16% for all three parton distribu- 
tion sets. It is noteworthy that unlike what obtains in cc 
production at this energy, both the direct and resolved 
components from EPS09 are the largest. In addition the 
resolved component from FGSIO is also larger than that 
of MSTW08. This is due to the increasingly important 
role of antishadowing at this energy. Thus the EPS09 
total cross section is 2.2% larger than that of MSTW08. 
Due to the stronger shadowing and smaller antishadow- 
ing in FGSIO, the total cross section is approximately 
1.1% smaller. 

\/^NN — 5.5 TeV the resolved components are re- 
spectively 21.5% for MSTW08, 22.4% for EPS09, and 
22.9% for FGSIO. While the influence of shadowing can 
be seen in the direct components, both EPS09 and FGSIO 
yield larger resolved components, again due to antishad- 
owing. Thus the MSTW08 total cross section is reduced 
by 0.6% and 4.4% by EPS09 and FGSIO respectively. 

In Fig. [5] we show the corresponding rapidity distribu- 
tions. Considering the left-hand panels (^s„„ = 2.76 



in the total distributions. 

These features are also present in the distributions at 
y/T^ =5.5 TeV. Here also the clear influence of shad- 
owing in the rapidity interval — 2 < y < 2 in the direct 
component is reduced by the antishadowing present in 
the resolved component in the same interval. Thus the 
total rapidity distribution exhibit reduced sensitivity to 
nuclear shadowing in this interval due to this destruc- 
tive interference. Nevertheless the influence of shadow- 
ing is still apparent especially in the rapidity window 
— l<y<l. Thus this interval presents the best sensi- 
tivity to shadowing effects in bb production. As at the 
lower energy, there is a slight manifestation of antishad- 
owing in the intervals — 4 < y < — 2 and 2 < y < 4. 

Let us now consider the magnitude of the change in 
cross section with energy. In Table |V] we present the 
ratio of the total cross section at ^/s^ — 5.5 TeV to 
that at ^s„„ — 2.76 TeV for cc, Rcc, and for bb, R^^i. 
As is readily apparent from the Table, both Rcc and Rj^i 
vary slightly with the choice of parton distribution set. 
Despite this variation, to a good approximation i?cc ~ 
1.7 and R^^i 2.4 for all three parton distribution sets. 
These values are larger, especially for i?;,^ than obtained 
for both pp and pPb collisions. 
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FIG. 8. (Color online) Rapidity distributions of bb photopro- 
duction in ultraperipheral PbPb collisions at the LHC. In (a), 
(c), and (e) we show the direct, resolved, and total distribu- 
tions respectively at ys„„ = 2.76 TeV while (b), (d), and (f) 
denote the equivalent distributions at ys„„ = 5.5 TeV. In 
all cases long-dashed (MSTW08), solid (EPS09), and dash- 
dotted (FGSIO) lines correspond to no shadowing, moderate 
shadowing, and strong shadowing respectively. 

TeV), it is apparent that the influence of shadowing 
present in the direct component in the rapidity interval 
— 1 < y < 1 is almost totally negated by the antishadow- 
ing dominant in the same interval in the resolved compo- 
nent, with the result that in this interval total rapidity 
distributions have very little sensitivity to nuclear effects. 
For — 3 < y < — 1 and 1 < y < 3 the slight manifestation 
of antishadowing seen in the direct component persists 



D. Theoretical errors and ratio of cross sections 

From a consideration of the general structure of the 
quantities (cross sections and rapidity distributions) con- 
sidered in this study, three major sources of theoretical 
errors can be readily identified: 

• accuracy of the relevant expressions for photon flux, 

• higher-order corrections, 

• uncertainties in parton distributions. 

Due to technical reasons we have not attempted to esti- 
mate quantitatively the uncertainties on calculated quan- 
tities in the present study. A pragmatic approach to fa- 
cilitate shadowing determination is to compare photopro- 
duction in proton-nucleus and nucleus-nucleus collisions, 
where many theoretical uncertainties and systematic er- 
rors cancel (see for instance [1, [l3|)- Here we explore 
in a rather simplistic and nonrigorous manner such an 
approach for shadowing effects in cc and bb photopro- 
duction. We are primarily interested in relative effects, 
and as such issues like exact normalization and flux par- 
allelism have been ignored. 

In Fig. El we show the ratio a^^P^ / {2Q%(jpP^) as a func- 
tion of rapidity for both cc and bb in their respective shad- 
owing regions. In order to have nearly identical Lorentz 
factor 7 the PbPb cross sections, a^^^^, have been cal- 
culated at y/s^ = 5.5 TeV and the pPb cross sections, 
o'P^'', at ^/s^ = 5 TeV. It should be noted that in our 
calculations equal rapidity corresponds to equal photon 
energy. 
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Let us compare the shadowing effects from these ra- 
tios at a specific rapidity, say at y = 0. Thus for cc 
the MSTW08 value is reduced by approximately 12.5% 
and 23.2% respectively by the shadowing in EPS09 and 
FGSIO. These reductions are precisely what one obtains 
from the rapidity distributions of cc in PbPb collisions at 
y/s„„ =5.5 TeV. This is also true for bb where the ra- 
tio reduction values 3.5% (EPS09) and 7.1% (FGSIO) are 
the same as from the rapidity distributions. It thus seems 
feasible, at least in this simplistic case, that relative shad- 
owing effects are transmitted without appreciable loss. 
Of course, in order to ascertain to what degree higher or- 



der effects are cancelled in these ratios, one should carry 
out the analogous calculations at leading order. Further 
work along this line is in progress. 



IV. CONCLUSIONS 

In the present study we have considered next-to- 
leading order (NLO) photoproduction of cc and bb in ul- 
traperipheral proton-proton (pp), proton- lead (pPb) and 
lead- lead (PbPb) collisions at LHC. In addition to its 
value in probing several aspects of QCD dynamics, pho- 
toproduction of heavy quarks can also aid in elucidat- 
ing and constraining some components of nuclear par- 
ton distributions. Although the parton distribution de- 
pendence is linear and different modifications are super- 
imposed due to the integration over parton momentum 
fraction x, both cross sections and rapidity distributions 
for cc in PbPb collisions manifest appreciable sensitivity 
to shadowing around midrapidity and a slight sensitivity 
to antishadowing at more forward and backward rapidi- 
ties. Thus cc photoproduction offers good constraining 
potential for shadowing, and a somewhat less potential 
for antishadowing. Despite the fact that photoproduc- 
tion of bb is less sensitive to nuclear modifications than 
cc, the influence of shadowing is evident around midra- 
pidity, and it thus offers some constraining ability for 
shadowing. Both cc and bb total photoproduction cross 
sections and rapidity distributions in pPb collisions show 
little sensitivity to nuclear modifications and are there- 
fore useful in shadowing determination via cross section 
ratios. In addition their resolved components are appre- 
ciable, especially for bb and thus it seems feasible that 
they could be of some use in constraining photon parton 
distributions. 

We thank C.A. Bertulani and M.J. Murray for helpful 
suggestions. 
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